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ABSTRACT 

The  LINAC  at  the  Naval  Postgraduate  School ,  Monterey,  was  used 
to  accelerate  electrons  to  energies  ranging  from  52  to  92  MeV  in 
order  to  study  the  energy  distributions  of  high  energy  electrons 

before  and  after  passing  through  layers  of  tin,  gadolinium,  and 

2 
lead.  The  thickness  of  these  materials  ranged  from  0.8  to  5.9  g/cm  . 

The  most  probable  energy  losses  agreed  with  the  theory  of  Blunck 

and  Westphal  for  all  materials  used,  while  distribution  half-widths 

2 
agreed  only  for  absorbers  of  thickness  less  than  3.0  g/cm  .  The 

thickness  at  which  theory  and  experiment  began  to  exhibit  a  notice- 
able discrepancy  was  found  to  be  dependent  on  the  atomic  number  of 
the  material . 

Where  comparison  was  possible,  results  of  this  experiment  generally 
agreed  with  the  findings  of  similar  works  concluded  previously. 
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I.   INTRODUCTION 

The  theoretical  energy  distribution  for  an  initially  monoenergetic 
electron  beam  which  has  passed  through  an  absorbing  material  was  de- 
termined by  Blunck  and  Westphal  [1].  Their  work  modified  previous 
theoretical  calculations  of  Landau  [2],  Eyges  [3],  Bethe  and  Heitler 
[4],  and  Blunck  and  Leisegang  [5]  by  adding  radiation  losses  to  the 
previous  predictions  of  ionization  losses  for  electrons  when  passing 
through  thin  absorbing  layers.  The  distribution  assumes  energy  losses 
small  compared  to  the  incident  beam  energy.  The  theoretical  treatment 
is  presented  in  Section  II. 

Several  measurements  of  energy  losses  by  high  energy  electrons 
passing  through  various  materials  have  been  performed  in  the  10  MeV 
to  150  MeV  range.  Of  note  are  the  works  of  Breuer  on  aluminium  [6]; 
Bumiller,  Buskirk,  Dyer,  and  Miller  on  aluminum  [7,8];  Goodwin  on 
copper  [9];  and  DeLeuil  and  Raynis  on  aluminum,  copper,  and  lead  [10]. 
Except  for  Breuer' s  work  on  the  Darmstadt  linear  accelerator,  the 
experiments  were  performed  on  the  LINAC  at  the  Naval  Postgraduate  School , 

These  experimental  results  are  in  general  agreement  with  theory  for  thin 

o 
absorbers  (<  2  g/cm  ),  but  discrepancies  with  theoretical  predictions 

are  reported  for  the  thicker  absorbers.  Such  discrepancies  are  expected 

for  the  thicker  absorbers  since  the  energy  loss  is  no  longer  small 

compared  to  the  initial  beam  energy. 

In  this  thesis  the  energy  loss  measurements  are  extended  to  include 

tin  (Z  =  50)  and  gadolinium  (Z  =  64).  These  metals  were  chosen  to 
provide  experimental  results  for  materials  with  atomic  numbers  between 
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copper  (Z  =  29)  and  lead   (Z  =  82).     It  was  thought  that  analysis  of 
these  materials  would  provide  more  conclusive  information  on  the 
effects  of  atomic  number  on  the  agreement  between  experimental   and 

theoretical   half-widths  as  reported  by  DeLeuil   and  Raynis  [10]. 

2 
Thicknesses  used  ranged  from  1.485  to  5.940  g/cm     for  tin  and  from 

2 
0.814  to  4.831   g/cm     for  gadolinium;  nominal   beam  energies  were  52, 

75,  and  92  MeV.  Thick  lead  absorbers  (2.825  and  4.236  g/cm2)  were 

investigated  as  an  extension  of  the  work  of  DeLeuil  and  Raynis,  and, 

2 
in  addition,  their  thickest  aluminum  absorber  (5.574  g/cm  )  was 

investigated  as  a  continuity  check  between  the  two  experiments. 

The  energy  loss  distributions  are  characterized  by  the  most  probable 
energy  loss  and  the  half-width.  The  half-width  is  the  full  width  of 
the  distribution  curve  at  half  maximum.  These  quantities,  obtained 
from  experiment  and  theory,  are  the  basis  for  comparison  between  ex- 
perimental and  theoretical  results. 
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II.  THEORETICAL  CONSIDERATIONS 

The  Blunck  and  Westphal  theory  of  the  distribution  for  the  energy 
loss  of  a  beam  of  monoenergetic  electrons  in  passing  through  a  layer 
of  absorbing  material  assumes  that  the  energy  loss  Q  is  small  compared 
to  the  initial  beam  energy,  E. .  Let  W(Q)dQ  be  the  probability  of 
energy  loss  between  Q  and  Q  +  dQ,  and  X  be  that  portion  of  the  loss 
Q  due  to  radiation.  Hence,  the  ionization  loss  is  Q  -  X.  Considering 
these  two  loss  processes,  the  probability  of  energy  loss  is 

Q 
W(Q)dQ  =  /  WT(Q-X)W<;(X)dXdQ  (1) 

o  l  3 

where  W,  and  IaL  are  the  energy  loss  distributions  for  ionization  and 
radiation  respectively. 

The  Landau  equation  [2],  as  modified  by  Blunck  and  Leisegang,  is 
used  for  the  energy  loss  distribution  due  to  ionization.  The  distri- 
bution is  expressed  as  a  function  of  Landau's  dimensionless  parameter 
A  as  follows: 

,2 


L 


h2  I     2 

b   +   yn 


dA  (2) 


C  y 

WT(Q)dQ  =  *(x)  dx  =  z :  rS-S ^    exp 

f  +  y, 

—       E. 

where  A  =  ^  +  In   4-  -  1.116  (3) 

aK        aK 

The  terms  used  in  equations  (2)  and  (3)  are  defined  as  follows: 

C  ,  y  s  and  x     are  constants  given  by  Blunck  and  Leisegang  [5]  and 

are  used  to  fit  Landau's  distribution  to  a  sum  of  gaussian  functions. 
R  is  the  absorber  thickness  in  cm. 
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The  quantity  "a"  is  a  function  of  the  atomic  number  Z,  the  atomic 
weight  A,  and  the  density  p  of  the  absorber;  and  b(=v/c)  of  the 
electrons: 


:  0.154  Z  P 


~  rf  r       MeV/cm. 
A  f 


(4) 


The  quantity  b  is  a  correction  to  Landau's  theory  given  by  Blunck 
and  Leisegang  [5]: 


b2  =  3^0  z     Vm_ln 


aR 


m 


2  E. 


.  yi  -  *>. 


(5) 


where  the  summation  is  over  the  ionization  potentials  of  the  atomic 
electrons,  and  N  is  the  number  of  electrons  with  ionization  potential 

m 

Q  is  the  average  energy  loss  due  to  ionization   (no  radiation)  for 

electrons  of  incident  energy  E.s  and  is  given  by  Sternheimer  [11,12, 

and  13]  as  follows: 


-    V 


B  +  0.43  +  In  Ei  -  /   -  C  -  a  (X-,  -  log10p/mc) 


m 


MeV  (6) 


wh 


ere  t  is  the  thickness  in  g/cm  ,  and  the  constants  A  ,  B,  C,  X,,  a( 


and  m  are  parameters  of  the  absorber  material.  These  parameters  for 
tin,  lead,  and  various  other  material  are  listed  in  reference  [13],  The 
parameters  for  gadolinium  are  not  listed  but  were  determined  by  the 
following  method. 

A  for  gadolinium  was  obtained  by  extrapolating  from  a  plot  of  A 
vs.  Z/A  for  the  absorber  materials  listed  in  reference  13.  The  result 
was  A  =  0.0624.  B  for  gadolinium  was  found  in  a  similar  manner  from  a 
B  vs.  I  (ionization  potential)  plot  where  I  was  determined  from  a  Z  vs. 
I  graph.  The  result  was  B  =  13.3.  C  was  determined  from  a  semi-log  plot 
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2 

of  I  A/pZ  vs.  -C,  and  the  result  was  C  =  -6.80.  The  values  for  a  , 

X, ,  and  m  were  more  difficult  to  determine,  but  fortunately  they 

contribute  little  in  determining  Q.  These  constants  were  numerically 

m 
determined.  The  values  for  the  term  a  (Xn-log-,0(p/mc))   were 

calculated  for  tin  (Z  =  50)  and  tungsten  (Z  =  74)  in  the  p/mc  range 
(p  is  momentum,  and  m  is  the  rest  mass  of  the  electron)  of  the  energies 
used  in  the  experiment.  Using  Z  as  a  basis  of  interpolation,  the  cor- 
responding values  for  gadolinium  (Z  =  64)  were  estimated  for  two 
representative  p/mc  values  (100  and  150).  The  value  for  X,  was 
arbitrarily  set  to  be  3.0,  and  the  a  and  m  values  were  determined 
from  a  simultaneous  algebraic  solution  based  on  the  interpolated  values, 

The  results  were  a„  =  0.418  and  rn  =  2.10.  The  error  in  the  determina- 

s  s 

tion  of  Q"  using  these  values  was  estimated  to  be  1%  or  less. 

For  Ws,  the  energy  loss  distribution  due  to  radiation,  Blunck  and 
Westphal  [1]  give: 

Ws(Q)dQ  =  BaR(Q/E.)aR  ^  (7) 

where 


a  =  1.4  x  10"3  ^ 


4/3  In  1&-   +  ]/9 


cm"1        (8) 


and  B  is  a  normalizing  factor  =  ,   R  +  , \  . 

The  distribution  of  total  energy  loss  according  to  Blunck  and 
Westphal  is  obtained  by  putting  equations  (7)  and  (2)  into  equation 
(1)  and  performing  the  required  integration.  The  result  is  the  energy 
loss  distribution  for  a  single  electron  of  incident  energy  E. .  For 
comparison  of  theoretical  and  experimental  values,  this  distribution 
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function  was  used,  with  corrections  to  account  for  the  finite  energy 
width  of  the  incident  electron  beam.  This  treatment  is  described  in 
Section  IV  and  Appendix  C. 
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III.  EXPERIMENTAL  PROCEDURE 

The  LINAC  of  the  Naval  Postgraduate  School  was  used  to  obtain 
electron  beam  energies  from  52  MeV  to  92  MeV.  This  beam  was  elastically 
scattered  at  90°  from  a  thin  (3.5  mil)  aluminum  scattering  foil,  passed 
through  the  absorber,  and  finally  analyzed  by  a  120°  magnetic  spectro- 
meter described  by  Kenaston,  Luke,  and  Sones  [14]. 

This  general  experimental  arrangement  was  similar  to  that  used  by 
Miller  [7,8]  and  DeLeuil  and  Raynis  [10]  with  the  exception  of  the 
removal  of  a  3.5  mil  aluminum  window  and  the  installation  of  a  ten 
channel  coincidence  counting  system.  These  changes  are  discussed  below. 

The  aluminum  window  removed  was  located  just  before  the  scattering 
chamber.  Since  the  presence  of  the  window  caused  a  broadening  of  the 
incident  electron  beam  distribution,  its  removal  decreased  the  energy 
width  of  the  incident  beam  for  this  experiment,  as  compared  to  previous 
works  done  at  the  NPS. 

The  ten  channel  counting  system  consists  of  ten  front  counters  and 
a  single  backing  counter  operated  in  coincidence  with  the  front  counters. 
The  entire  ten  channel  system  has  an  energy  spread  of  about  3%. 

The  absorbers  were  positioned  approximately  3  cm.  from  the  scattering 

foil  as  recommended  by  DeLeuil  and  Raynis  [10].  Absorber  thicknesses  for 

2 
tin  were  1.485,  2.970,  4.455,  and  5.940  g/cm  ;  thicknesses  for  gadolinium 

2 
were  0.814,  1.610,  3.221,  4.026,  and  4.831  g/cm  ;  absorber  thicknesses 

for  lead  were  2.825  and  4.236  g/cm2. 

A  thick  aluminum  absorber  (5.574  g/cm  )  was  used  to  correlate  the 

results  of  this  experiment  to  the  results  of  DeLeuil  and  Raynis  [10], 


17 


which  were  obtained  before  the  ten  channel  countin.,  system  Wa  .  called. 
The  results  agreed  within  experimental  accuracy. 

Since  the  electron  beam  scattered  into  the  absorbing  matei    was 
not  monoenergetic,  the  energy  distribution  was  determined  both  before 
and  after  passing  through  the  absorber.  Various  thicknesses  of  absorber 
material,  including  no  absorber,  were  positioned  on  a  remotely  controlled 
ladder  device.  For  measurement  of  the  distribution  before  passing  through 
the  absorbers  (called  a  zero  peak),  the  ladder  was  positioned  such  that 
the  beam  passed  only  through  the  scattering  foil.  After  measuring  the 
zero  peak  spectrum,  the  various  absorbers  were  moved  into  the  beam  to 
determine  the  energy  distribution  of  the  beam  after  passing  through  each 
absorber. 

In  the  case  of  tin  and  lead  all  absorber  thicknesses  could  be  measured 
without  turning  oft  the  accelerator  and  thus  possibly  altering  the  char- 
acter of  the  beam.  This  was  not  possible  for  gadolinium.  Because  of  its 
cost,  only  a  small  quantity  was  purchased.  Hence,  only  two  gadolinium 
absorbers  could  be  run  without  turning  off  the  accelerator  and  rearrang- 
ing the  gadolinium  on  the  ladder.  Since  the  beam  character  could  change, 
a  zero  peak  measurement  was  made  whenever  absorbers  were  replaced. 

The  data  represented  the  number  of  electrons  detected  by  the  coin- 
cidence counting  system  at  the  exit  of  the  magnetic  spectrometer.  A  down- 
stream Secondary  Emission  Monitor  was  used  as  a  standard  for  normalization 
purposes,  in  that  each  data  point  corresponded  to  a  given  integration 
current,  i.e.,  a  certain  number  of  electrons  passing  through  the  absorb- 
ing material . 
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IV.  TREATMENT  OF  DATA 

The  data  reduction  for  the  ten  channel  counting  system  is  standard 
for  multi -channel  systems  and  is  on  file  in  the  NPS  LINAC  computer 
library.  In  this  reduction  process,  three  characteristics  of  the 
counting  system  are  required.  The  characteristics  are: 

(1)  the  energy  spread  of  the  counting  system, 

(2)  the  energy  seen  by  each  front  counter, 

(3)  the  relative  efficiencies  of  the  front  counters. 

With  this  information  the  actual  energy  corresponding  to  the  front  counter 
data  is  calculated.  Counting  rate  and  background  corrections  are  also 
performed  during  this  computerized  process. 

In  the  lower  energy  regions  of  the  spectrum,  where  all  counter  read- 
ings should  be  nearly  equal,  a  counting  discrepancy  was  noted  in  that  the 
lower  energy  front  counters  read  consistently  higher  than  those  on  the 
higher  energy  portion  of  the  counting  spectrum.  This  effect  occurred  at 
energies  lower  than  required  to  obtain  distribution  half-widths  and  thus 
did  not  affect  the  data  reported.  However,  it  made  confirmation  of  the 
shape  of  the  distributions  at  \/ery   low  energy  (<  20  MeV)  impossible. 

The  electron  beam  energy  is  not  monoenergetic  as  required  by  Blunck 
and  Westphal  theory.  To  compare  experimental  and  theoretical  results, 
the  zero  peak  data  must  be  used  to  modify  the  theoretical  predictions. 
With  no  data  treatment,  the  theory  will  predict  half-widths  which  are 
too  small  and,  in  the  case  of  asymmetric  zero  peaks,  erroneous  most 
probable  energy  losses.  Removal  of  the  aluminum  window  noted  in  Section 
III  narrowed  the  zero  peak  to  the  extent  that  it  could  be  treated  as  a 


19 


symmetrical  gaussian  distribution.  The  IBM  360  computer  of  the  MPS 
was  used  to  unfold  the  zero  peak  energy  distribution  into  the  theory. 

A  histogram  method  described  in  Appendix  C  was  used  for  this  unfolding. 
The  computer  program  used  to  accomplish  the  unfolding  is  Appendix  D. 

In  the  comparison  of  experiment  with  theory,  the  measured  incident 
energy  distribution  was  unfolded  as  previously  described,  properly 
normalized  to  the  experimental  data,  and  plotted.  From  these  plots 
the  theoretical  half-widths  and  most  probable  energy  losses  were 
determined.  The  experimental  data  were  plotted  along  with  the  correspond- 
ing theoretical  curve,  and  measurable  parameters  were  compared.  The  curves 
are  shown  in  figures  1  through  38  of  Appendix  B. 
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V.  RESULTS  AND  OBSERVATIONS  * 

The  values  predicted  by  theory  and  the  experimental  results  for  the 
most  probable  energy  loss  and  half-width  for  tin,  gadolinium,  and  lead 
are  shown  in  Tables  I,  II,  and  III  respectively.  In  these  tables,  BW 
refers  to  Blunck  and  Westphal  theoretical  predictions,  Q  represents 
the  most  probable  energy  loss,  and  Hw*  is  the  half-width.  Data  are  given 
for  predictions  with  and  without  beam  folding. 

On  comparing  the  theoretical  predictions  with  beam  folding  and  the 

experimental  results,  it  is  seen  that  good  agreement  exists  for  target 

2 

thicknesses  less  than  about  3.0  g/cm  for  both  half-width  and  most 

probable  energy  loss.  At  thicknesses  greater  than  this,  the  agreement 
in  most  probable  energy  loss  is  still  reasonably  good,  but  there  is  poor 
agreement  in  half-width.  The  greatest  variation  in  most  probable  energy 
loss  is  8%  while  the  average  variation  is  3%.     These  percentages,  while 
averaged  over  all  material,  are  typical  and  not  material  dependent.  The 
average  half-width  variation  is  9%  where  half -widths  were  obtained.  As 
can  be  seen  from  Appendix  B,  figures  2,4,  and  6,  the  theory  does  not 
predict  a  half-width  for  a  certain  target  thickness  for  each  element. 

This  cutoff  thickness  is  smaller  as  the  atomic  number  increases.  For 

o 
example,  the  theory  would  predict  a  half-width  for  5.6  g/cm  aluminum, 

2 
but  not  for  4.2  g/cm  lead. 

It  is  concluded  from  these  results  that  the  Blunck  and  Westphal  theory 

predicts  correctly  the  most  probable  energy  loss  for  all  targets  used  and 

predicts  satisfactorily  the  half-widths  for  the  thinner  targets.  Thin 
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targets  can  be  defined  empirically  as  those  satisfying  the  relation 

1  /3 
T(Z)    <  13.5,  where  T  is  the  absorber  thickness  (g/cm2)  and  Z  is 

the  absorber  atomic  number. 
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APPENDIX  C  -  BEAM  FOLDING  TECHNIQUE 

The  theory  of  Blunck  and  Westphal  [1]  assumes  that  the  beam  of 
electrons  striking  the  absorber  is  monoenergetic.  The  beam  of  electrons 
produced  by  the  NPS  LINAC,  or  any  other  linear  accelerator,  for  that 
matter,  has  a  finite  energy  spread  about  a  maximum  or  most  probable 
energy  point.  The  fact  that  monoenergetic  electrons  are  not  avail- 
able to  strike  the  absorber  must  be  taken  into  account  in  computing 
theoretical  predictions  if  a  meaningful  comparison  with  experimental 
results  is  to  be  made.  This  has  been  done  here  by  a  technique  termed 
"beam  folding". 

Beam  folding  is  accomplished  by  a  number  of  well  defined  steps. 

Energy  distribution  curves  may  be  approximated  by  histograms.  These 

histograms  consist  of  a  series  of  "bins"  of  area  W(Q)aQ,  where  Q  is 

the  energy  loss  and  W(Q)aQ  is  the  probability  of  loss  between  Q  and 

Q  +  aQ.  Thus,  each  bin  has  an  "address",  Q,  on  an  energy  coordinate 

scale.  To  accomplish  beam  folding,  the  beam  distribution  must  be  known. 

This  is  observed  experimentally  and  approximated  in  the  computer  by  a 

gaussian  curve  of  appropriate  half-width.  The  energy  at  which  the 

maximum  occurs  is  established  by  the  energy  of  the  beam  and  the  magnitude 

of  energy  loss  incurred  by  elastic  scattering  of  the  beam  as  it  impinges 

on  the  thin  aluminum  scattering  foil  prior  to  striking  the  absorber. 

The  width  aQ  is  then  selected  for  the  predicted  distribution.  This  must 

be  a  small,  but  finite  number  where  numerical  techniques  are  to  be  used. 

This  same  width  is  used  to  break  the  beam  distribution  into  a  histogram. 

This  is  illustrated  in  figure  1(a)  of  this  appendix.  The  reason  the 

same  width  is  used  for  both  distributions  is  to  facilitate  computer 

programming. 
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(a)  Beam  Distribution 


(b)  Energy  Distribution  of  E 


(c)  Energy  Distribution  of  E,- 


(d)  Total  energy  Distribution 


Energy 


Energy 


Energy 


Energy 


Figure  1 
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Each  bin  of  the  beam  distribution  is  now  treated  as  a  monoenergetic 
beam  with  energy  commensurate  with  its  center  location  on  the  energy 
scale.  Each  bin  also  has  a  definite  magnitude,  or  weight,  with  the 
magnitude  of  the  center  bin  being  unity.  The  formulae  of  Blunck  and 
Westphal  is  now  applied  to  each  of  these  beams  and  an  absorber  distri- 
bution curve  results  for  each,  with  a  maximum  amplitude  proportional 
to  the  height  of  the  appropriate  beam  distribution  histogram.  Each  of 
the  absorber  distribution  curves  thus  obtained  may  be  thought  of  as 
being  plotted  and  added  to  previously  determined  curves,  using  the 
energy  of  the  electrons  as  a  correlation  point,  as  depicted  in  figures 
1  (b),  (c),  and  (d)  of  this  appendix.  This  is  accomplished  on  the 
computer  by  adding  the  contents  of  each  bin  of  the  same  address  and 
plotting  the  cumulative  total.  The  bin  of  address  E  is  depicted  in 
figures  1  (b),  (c),  and  (d).  Note  that  the  beam  distribution  for  histo- 
gram E4  is  centered  over  the  maximum  point  of  the  distribution.  This 
is  important  as  a  false  picture  could  easily  be  presented  if  the  histo- 
gram were  not  symmetrical  as  the  distribution  would  then  appear  skewed. 

Since  the  half-width  for  a  particular  absorber  increases  as  the 
thickness  (g/cm  )  increases,  the  beam  folding  technique  is  affected  in 
that  thick  target  energy  distribution,  whose  half-width  is  large,  will 
be  affected  much  less  by  the  finite  beam  distribution  than  will  the 
energy  distribution  of  a  thin  target  whose  predicted  half-width  is  not 
more  than  2  or  3  times  the  size  of  the  beam  distribution  half-width. 
Thus,  in  selecting  bin  width  (aQ),  it  is  desirable  to  make  the  width 
proportional  to  the  target  thickness  so  as  to  apply  the  beam  folding 
technique  in  greater  detail  to  thin  as  opposed  to  thick  targets.  The 
computer  program  (Appendix  D)  applies  the  target  thickness  in  direct 
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proportion  to  establish  a  bin  width  for  subsequent  beam  folding 
calculations.  Thus  aQ  is  much  smaller  for  the  thinner  absorbers 
compared  to  thick  ones. 
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APPENDIX  D.     Computer  Program 


C    PURPOSE:    THIS    program    ACCEPTS    p<\damptF"<;    F^3m    cMrPr,v    |  "> «; c; 
C        fxpcrt  mpnts   roNnurTFn   on   the    \  injaf    a  mo   wtli    compute 

f.    THEORETICAL  HALFWIOTfi  AMD  MOST  PRORABLF  =NCRGY  LnSS  AVH 

C    WILL  PLOT  ROTH  THE  FXP  ER  I MFNT  AL  AND  THEOR  FT  I C  AL  cNr°r,v 

C    DISTRIBUTION  FOR  THF  EXPFrimfnt, 

C 

C  SCOpF:  THE  PROGRAM  TS  SPT  UP  TO  MAKP  C  OMPUT  A  T I  on  <;  prL4T]Vn 

C    TO  CXPFR  TM<=MTATin\!  WITH  ALUM*  NUM  ,  C  oppcr  t  t  I  N,  GADO!  !  m.  i  ||m 

C    AND  LEAD, 

C 

C  METHOD:  THE  PROGRAM  PLnT^  A  Theory  HPVE  RASED  nN  THE  CTT 

C    OF  GAUSSIAN  CURVES  PPRFORMcn  Ry  RHJNCK  AND  WPSTPMAL  rOR 

C    THF  PREDICTED  i=NcRGY  0  I  S  T°  T  BUT  I  ON  CURVES,   THF  TNTEGRA- 

C    TION  n\/FR  PVj^PGY  LOSS  IMDTC.ATFD  IN  THF  4  PHR  c  yc  kit  ]  hnc  n 

C    PAPER  IS  pFRFORMPD  ON  TMP  CQMDUTPR  USTNO  "*  ?  POINT  ^.  /\l  I S  S 

C    OUAOPATUPP,   TWO  GRAPHS  ARF  PL°TTFD;  ONr  "FOQ^^FS  THF 

C    R-W  PREDICTION  WITH  A  MONFN=R  Gc  T  T  C     rc^m  oc  FLpCTRONS  AMI 

C    THF  OTHPR,  WHICH  ALSO  "AS  rx  P  FP  I  mpnjT  A!   PF^IILTS  THcprnM, 

C    RFPR^SPNTS  THE  SAMF  THcr»RY  WITH  A  POLYFMC  R  GE  T I C  RnAM 

C    DISTRIRUTION  FOLDFO  IN, 

C 

C  THIS  PROGRAM  IS  WRITTEN  USING  FORTRAN  TV  IANGUAGC  ANO  tmc 

C    NPS  COMPUTFR  FACILITY  PLOTTING  PACKAGE,   TOTAI   RUNNING 

C    TIMF  IS  ABOUT  25/T  "MINUTES,  WHFR  F  T  IS  THF  ARSDRBrR 

C    THICKNESS  IN  OM/SO  CM, 

C 

IMPLICIT  REAL*P( A-H.O-Z) 
C 

C  FSTABLISH  TITLES  FOP  PLOTS* 
C 

REAL*R  Tl(4)/«  ENERGY  DISTRIBUTION  •/ 
REALMS  T3(4)/«  ALUMINUM  ABSORBER  •/ 
RCAL*8  TM4)/»T:  GM/CM  ej:  m^V'/ 
RFAL*B  T5(M/«  UNFOLDFD  THFORY  •/ 
REAL*8  T6(4)/«THF0RY:  FXPErtmcnta|  DATA:X»/ 
RFAL*8  T7(4)/«  LEAD   ABSORBER  •/ 

REAL*8  TB(M/»  TIN  ABSORBER  •/ 

REAL*P  TQ(4)/«        GADOLINIUM   ABSORBER        t/ 
PFAL*8  TlOC^t/1         COPPFR   APSHRRFR         t/ 
D  IMF  MS  I  ON  «=(«50)  ,  WF  (  5  0)  ,WQT(  ?O'MtHW(30'M,XX(?0'M, 
lWUFnO^)fXUF(^n^),ExyCOO)trxY(POO)tiyY(?CO)fT'>(^) 
EXTFRNAL  FCT,RFAMt WIOMtWRAD 
D01200  J=1,eo 

E( J)=0 

1?00  WF(J)=C 

D01201  J=1,?oq 

FXX( J)^0 

FXY( J)=0 
1201  IXY(J)=0 

DO1202  J=I t^OO 

WOT( J)=0 

HW( J)=o 

XX( J)=0 

WUF( J)=0 
120?  XUF(J)=o 

C 

C  PFAO  TN  PAPAMFTE^S  FOR  EXPERIMENT-. 

C  FB  BEAM  cnjproy 

C  7  ATOMir  numrer 

C  A  ATOMIC  WEIGHT 

r  T  TARGFT  THK^NESS:  C,m/sq  CM 

C  07EP0  SPFfTPnMFTPP  INDICATION  OF  LOCATION  np  rc^M 

r  DISTRIBUTION  (Mcy) 

C  CTSNnR   NUMBER  TO  WHICH  eXppdjmfnTAL  COUNTS  WILL  RE 

C  N0RMALI7EH 

C  HT  HALFWIDTU  OF  BEAM  DISTRIBUTION 
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RCAP(5, 11 )     EB, Z,A, T,OZERO,CTSMOR,HT 
C 

C     PTVIDE     7CRP     PEAK     INTO    "RIMS"     ^    WIPTH    DE|_X,        BT*1     WIPTH 
C  INCREASES    WITH    TAPGFT    thtCKM«=SS, 

C 

DCLX=.03J'T 

IF(T,LT0!„  )    PELX=,  Of>*T 
C 

C  OAVF  IS  THF  RFCPTL  LOSS  OF  ^t crTP1MS  STRIKING  AllJMINUM 
C    NUCIETT  OF  TH«=  SCATTER tng  fpil,   IONI7ATIPN  LOSS  TN 
C    SCATTFRF"  IS  NFGlECTEO„ 
C 

QAVE=FB*EB/(931,4*A) 

c 

C  rnRR  IS  thp  rnPRrr.TinM  to  be  aopco  to-  exper  I  mcnjtal  "=mci?ry 

C    TO  OORRFCT  FOR  SPrr.TRHMFTPR  CRPPR, 

c 

CPRP=0?FPP-QAV,= 
12  FORMAT  I  5(1  X,Pfl.-,?,l  X,  16)  ) 

c 

C    READ     IN    UP     TP    ?^P     FX»FR  TMFMTA!      POINTS;     CXX     IS    ENERGY     («rVI 
C  AND     IXY    TS    COUNTS,        THERE    MUST    RE    ^0    OATA    CAROS    AMP    the 

C  LAST    E XX    MUST    BE     7Fon, 

C 

READ (5,1?)     (EXX(J)  , TXY(J), J= 1,200) 
C 

C    ESTABLISH    PAPAMFTERS    FPP    SCATTFRINO    FOIL-, 
C 

FOTL=,0">4 

ZFOIL=13„ 

KSIG=0 

00    1     J=l,200 
1     EXY( J)=TXY( J) 

M  =  0 
C 

C  pI  IS  VALUE  OF  BEAM  ENERGY  INCIDENT  PM  TARGCT, 
C 

FI=EB-PAVF 

TEMP=PI 


C 

c 
c 
c 
c 
c 
c 


THF  FOLLOWING  portion  pF  THE  PR3GRAM  UNFOLDS  BFAM  OtsTdT- 
BUTTON,   ICENT  IS  AM  TNTE0E3  VALUE  ASSIGNED  to  tHc 
HISTOGRAM  CFNTerfo  nyr-o.    the  PFAK  OF  THE  8cAm  DTSToip.ii- 
TIOM,   FACH  HISTOGRAM  IS  DCMOTPO  RY  ITS  CMCRGY  e  AMP 
MAGNITUDE  RELATIVE  Tn  th^  p^AK  OF  THF  n  I STR  I  BUTT  PM,  WF„ 

ICPNT=(HT+PELX)/PELX 

X  =  EI-ir.FNTJPELX 
20?    X=X+PFLX 

IF(X.GT„ (EI+HT) )  GPTO  1003 

M  =  M+1 

F(M)=X 
C 

C  FUNCTION  BrAM  REPLICATES  THC  P T STo TRUT I PN  or  BEAM  cvjcp^y 
C  AS  A  SIMOLF  GAHSSTAM  CURVE  BASEO  PM  THE  EX  Pc  R  T  ^rMT  ALl  Y 
C  PBSFPVFO  HALFWIOTH  ANO  TS  IPCATCP  RELATIVE  Tn  THE  r^am 
C    ENERGY  BY  THE  COMPUTED  AVERAGE  RECOIL  LOSS, 

C 

WF(M)  =  REAM( EI, HT,X  ) 

GOTP?0? 
1003    WRITER,??)     0AVE,HT,EPH  t7EPjL,EI 

DO    100f    L-]  ,M 

DX  =  0 

N  =  0 
10C^    N  =  N+1 
C 

C  USING  EACH  HISTOGRAM  AS  \     "BCAM",  THF  EP|_LOWTNO  PnuTl^j  OF 
f    the  PRPG^AM  fpMPUTES  AM  FNFRGY  niSTRTBUTIOM  F^R  FAfH 
C    OF  THESE  "RF\ms«  AMD  STORES  THFM  IN  "BINS"  LARFLEP  RV 
C    EWFPGY,  XX,   The  CUMULATIVE  AMPUNIT  PE  F!_eftppm  COUNTS  IM 
C    THE  BIN  TS  WQT,   0  TS  THE  ENERGY  LPSS  PVER  WHICH  THr  R-W 
C    PPFDTCTIPN  IS  INTEGRATED  AMD  OX  IS  THE  AMOUNT  BY  WHICH 
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c 

r 
C 


C 

r 

C 

r 

C 

c 


0  CHANGES  FOP  EACH  Sljrr.FSSTV"  TN 
TS  CAPPED  OUT  IN  SUBROUTINE  DOG 

O=3«0*TM7/13,  )'  •>,  ^-r>X 
IP(Q9LT^EI  JGOTO1O04 
N  =  0 

GDT01005 
100^  IF(0,LTon)r,OTni0O6 
EI=F(L) 

VALIM=„?5*Q 

CALL  D003?(0, t XL    , 

CALL  nOG3? ( XL, VALI M, FfT, 

CALL  DOG^?( VALIM,0  ,f=CT, 

WQ=ZA+ZB+7C 

WQTCL+N-1 )=WQT (I  +  N-1 )+WO*WF(L) 

XXU +  N-1 )=F(L)-0 

J=L+N-1 


TFT,o  ATION. 
3?* 


TMTCGP  ATI™) 


FCT,7A,F 

7B,P 
ZCtE 


OtTt  A,Z> 
Of Tf A, 71 
0»T, A, 7 ) 


THE  DISTRIBUTION  fna 
CFNTpR  OP  THf  BFA^ 
WUF  A\in  XI  IF  AS  THE 
RESPECTIVELY, 


TH<=    HISTOGRAM 
niSTP  TBIJT  irjw 
UNFOLDED    CDUN 


nnpR  PSPHMDTNG 
IS  PPMoy/FD  A  NO 
TS    AND    ENERGY 


Tn  thp 
ST1Rr0  IN 


1006 


C 

c 
c 

c 
c 


IF(!  ,NF„  ICFNT)G0T01005 

WMF( M) =W0 

XUF<N)=FI-0 

K  =  N 

GOTH  1005 

CONTINUE 

On  1717  1=1,7 


FLIMINATP  LOW  NUMBEBFD  "RTNS"  AS  T 
LEO  ANO  ARF  NOT  P E PP E SPNTATI VF  1 
DISTRIBUTION, 


H<=y  ARF  ONLY 
F  THE  CUMULA' 


PARTLY 
I  Vp 


F!!.- 


C 
C. 

c 

c 
c 


WUF(  I)=W!)F(8) 
1717  WOT( I)=WQT( «) 
Dn  6  1=1,4 
6  T3(  I  )  =  T6(  I  ) 

Goth  77P 
777  D0778  N  =  1,K 


PPHCPSS  THP  UNFOLDED  THEORY  CURVE. 
TO  ZEPO  ALLOWS  PLOTTING  HP  THIS 
EXPERIMENTAL  DATA„ 

WQTf N)=WUF<N) 

X  X  (  N  )  =  X ! J  F  (  N  ) 
KSIG=1 
J  =  K 
778  CONTINUE 
EXX( 1)=0 

nn  5  i  =  i  ,4 

5  T3( I >  =  T5t I  ) 

77Q  WQTMX=0 

001007  N=3,J 

DFTrRMT\]P    7Hn     MAXIMUM    pniNT    nn     th^ 
AND    THUS     cSTAPLISi-<    THC    ^OST    PPOB 


SETTING    EXXC  M     cOUAL 
TH^OPY    CUPVE     WITHOUT     ANY 


HI STP I  RUT  ION,     WOT  MX, 
ABLE     >=NcROY    LOSS,     QP, 


TF(WOT{ M)„  1  cnwOT(N-l  )o  ANO, WQT(M-i  )  , GE . WOT ( N-? )  ) 
1GOT01003 
GOTn   loo*7 

10O8     IF(WOT(  M-l  )  ,|  T,  W0TMX)r,OTP    10^7 

WQTMX=WQT( N-l ) 

EI=TEM° 

QP=F  T-XY(N-T  ) 
1007  CONTINUE 

i  NORVALIZF  ALL  POINTS  OP  DTSTPIBUTTnN  AND  ST^PF  TN  HW, 
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1 


c 
c 
c 

c 
c 
c 

c 
c 
c 
c 
c 

r 
C 


on  ioqc  m=i ? j 

HW( M )=WOT( M) 7WQTMX 
^C°  CONTINUF 
J  =  J-1 

HWI=0 
DFTFRMIMP  HALFWIDTH  hf  THc  TH^HRPtICM  n  I  STR  T  RUT  I  ni\)  t  HWA , 

no  1010  Mai, j 

IF(HW(  M*i  )  ,GT.  ,*,  ANP,  HW(  Ml  .,  I  T,  > 5  )HWI =XX ( M ) +(  ,  5- 
1HW(M)  )  '■•nr?LX/(HW(  v+1  )-HW(  V)  ) 
IF(M,FQ,  J„  4MD,HWI„  Pf),0  )r,nTf)iiq 

GnTni2l 

1 \o     WPITP(6,?0) 
HWA  =  n 
GOTOIOi  n 
?0  FORMAT  (/  ,?X,  •  NO  HALrWIDTH  OBTAINFrv) 
1  ?  1  I  c  (  HW  (  M  )  o  G  T  ,  „  5  ,  A  MO „  HW  (  m  +  ]   )3LT>0«5)HWA  =  XV{M)-f(HW{M) 

1-*5)*DELX/<HW(M)-HW(M+] ))-HWI 
010  CONTINUE 

ESTABLISH  CUTOFF  SIGNAL  COR  PLOTTING  SUBROUTINE-! 

XX( J+l )=0 

ESTABLISH  CORRECT  GRAPH  TITLES  FPR  ABS^RFR  USED, 

IF(Z.GT,13, l)G0Tn2Q 
DO  3  1=1 ,4 
3  T2(I  )  =  T3U  ) 
GOT01011 

?P  IF(Z„GT,?Q^! )GHTH50 
DO  101=1,4 

10  T?( I  )  =  Ti 0(  I  ) 

GOT01011 

50  IFC7»GT,50.3l)GOTn64 

On  8  1=1,4 

8  T2<I )=T9  ( I ) 

GOTH in ii 

64  IE(Z«GT*A4,nG0Tn8? 
00  9    1=1,4 

9  T2(I )=TQ  ( I ) 
G0T0101! 

p?  no  7  1=1,4 

7  T2< I )=T7  ( I ) 

SEND  DATA  TO  PLOTTING  SUBROUTINE  FOR  PROCESSING.. 

Oil  CALL  GRAPH(XX,WQT, EXX,EXY,T4,T3,T2,T1,HWA,EI,T, 
lCOPR,CTSNHRtOP) 

PRINT  RESULTS  OF  COMPUTATIONS  IN  TABULAR  F0RMo 

WRITE(6,22)  0P,HWAtT,7t'rI 
IFfKSIO^O^)  GOTO  ^77 

11  FORMAT  (7^1.0,0) 
22  FORMAT{ /,2X, •  0°:  •,F7,4,» 

1'  T(GM/SO-CM):  »,F6-,4,«  7: 
2,F6«2,/) 

STOP 

END 


HALC  WIDTH:  •  ,F*.  i 


1  INCIDENT  EN-PGY:  • 
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c 

C  PI) 

c 
c 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


SUBPntJT  INF  nOG32(XLfXUtFCT,Y,F  I,0,T,Dt7) 

RPOSF:  THIS  SUBROUTINE  pfpftpmS  THF  INTEGRATION  RFQUIPFO 
IN  THF  R-W  THEORY. 

IMPL 
A  =  ,* 
B  =  XU 
C  =  «4 
Y=,3 
KFCT 
C  =  *4 

Y  =  Y  + 
!  (FCT 

C  =  ,4 

Y=Y  + 
1(FCT 

C  =  *4 

Y=Y  + 
1  (FCT 

c=,,4 

Y  =  Y  + 
1  (FCT 

C  =  *4 

Y  =  Y  + 
1  (FCT 

C  =  ^ 

Y  =  Y  + 
1  (FCT 

C=.3 

Y  =  Y  + 
1  (FCT 

C  =  ,3 

Y  =  Y  + 
1  (FCT 

C  =  ,2 

Y  =  Y  + 
1  (FCT 

C  =  ,2 

Y  =  Y  + 
1  (FCT 

C  =  ,2 

Y=Y  + 

KFCT 

C=.l 

Y  =  Y  + 
1  (FCT 

C  =  „l 

Y  =  Y  + 
1  (FCT 

C  =  .7 

Y  =  Y  + 
1  (FCT 

C  =  ,2 
Y=R* 
KFCT 
RrTU 
END 


ICTT    P 

raL*-8(  A-H,n-7  ) 

oom xu 

+  XL  ) 

-XL 

Op  A3"»  Q 

3C9?474O7Pn0+  R 

50Q7oe;liA4^c;oiQ5n-')* 

(( A+0) 

,7,n,T,FI,0)+FCri( 

P?«05T 

^57  7?f  34]  700 ''R 

•8137! 

Q"M  A646?  P3RO-?  <r 

((A+C) 

♦  7,DtT,c  Tf 0)+FCT  (  ( 

823811 

?77Q77^3?  ?oo*R 

ol  ?6C"S 

n*?  654  6"*  1^300-1  * 

((A+C) 

iZtOtT.EIf 0>+FCT( 1 

^,7^5  30^7Q6PP^Q«  4HO  *R 

«!71?6 

Q^l  4^*^!  07170-1  * 

(( A+C) 

,7.  0,TtFTrO)+FC,T(  ( 

4R16^577RR"*0?69An(y~B 

,21417940  011 H  3340D-1* 

((A+C) 

t 7f n,T,FI, 0)+FCT( ( 

246R8R 

'V  366?  R49900V'R 

„?K4Q9 

n?q5^i  i  flpnsgn-]  •- 

((A+C) 

»  7,0,T,FI tQ)+FCT  (( 

97?  M  R 

0795-3071  ?onn  -n, 

,?Q34?0£*73Q?f,77  74n-'i*- 

(  (  A  +  C  ) 

,7.n,TtCTt0)+FCT ( ( 

66^C»10S°37<"-1  44R40O*R 

-,^7c5]11113Rc»lR^r)2^n-l^ 

((A+C) 

t7,0,Ttcit0)+Fr.T  (  ( 

315221 

334651  0760m*R 

,-5(c,17?pQ7  0S4^?^?53r»-l»- 

(  (  A  +  C  ) 

.7,D»T,EI,0)+FCT ( ( 

93P5787»620^116nr»*R 

0 3 90 96 

947893535 15^0-1  * 

((A+C) 

,7»ntTtFIt0)+FCT ( ( 

5344995^4661  14'7C0^,*?R 

.4165596?ll?47  3*7*n-l* 

((A+C) 

, Z,D,T,ET,0)+FCT ( ( 

106756 

3B06531  76700  ?;B 

a43R'>6  94£  qn??r^i9n^n-i* 

((A+C) 

♦7t0,TtFT,0)+FCT{ ( 

65°^4^ 

01141 ^6"*R?nO*R 

« 45 5 86 

O3qr>A7RRlO4^0-l4- 

((A+C) 

•7tDtT.FI. 0)+PCT(f 

1Q64  36R112  6^6R5400&R 

e46°??lQO5^r>4n??fl30-1* 

(( A+r) 

,7,0tT»EI,0)+FC^( ( 

723  598  07913982  50-1 ^R 

,47  RIO 

■2  AO0?°6'*7  4'^0-!  * 

((A+C) 

,7»6|TfF  I,  6)  +FC.T(  ( 

41  5  3  33 

?P4^  R69]  50Q-''  *  R 

(Y  +  „4R?79  0  44?67T6^QAnr)-l 

((A+C) 

, 7tOtT,Ei,o)+FCT( ( 

RN 

A-C) tz ,nf t,fi tQ) 


A-C) ,Z ,ntTtFI,0) 


A-C)t7,OtT,ET,0) 


A-C)  ,7  ,n,T,nfQ| 


A-C) ,Z tO» T,c!f0) 


A-C)  ,Zt0,T,FI,0) 


A-C) , 7,0, T, CT,Q) 


A-C)  ,Z tO,TfFI,0) 


A-C) tZ»n,T,Fi,o) 


A-C) ,7fD,T, cTf q) 


A-C)  »Z ,Df T, Fit  0) 


A-C) f7,0,T,EI,Q) 


A-C)  ,Z  ,OtTt<=It0) 


A-C)  ,7  fO,T,FI, 0) 


A-C)  »Z  ,0»T,fi, Q) 


A-C)  ,7tOtTtFI,0)  ) ) 
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REAL  FUNCTION  WI ON *B ( X , r I , T f 7 , A , 0  ) 

PURPOSE:  THTS  FUNCTION  CALCULATES  THAT  PART  nE  THE  enfrqy 
DISTRIBUTION  OF  AM  H  cf.T^N  OUT  Tn  L^^S  HF  emcqry  «Y 
IHMIZATI^N  WHILF  PASSING  THROUGH  AN  ABSORBER, 

IMPLICIT  R^AL*R(  A-H,0-Z) 

IL     ARC    THE     AVERAGE     I  ON  T  7  AT  T  nri    POTENTIALS     PEP     ^LTTRHM    FOP 
LEAD    FOR    EACH    SHELL,     BEGINNING    WITH    THE    K     ^HCLL, 

NPOTL     IS    Thf    NUMBER    °F     inNI7ATin\j    POTENTIALS    USED    r"iP    !c/\o 

TN    THIS    CONTEXT,     L=LEAO,    G=GADOL IN IUM,     S=TIN,     C=COPPER, 
AND    A=ALUMINUM-> 

REAL**  IL(  6)/.  0380 00, »  01 4700,,  003200, ->  000490, 
1*000076,.,  00  0002/ 

RFAL*R  NL(6)/2-»,8,,lRn,*?,,lR,,4,/ 

PCAL*R  IG(5)/,OSO?4t,no7=?5,or>^lQ9f,Or)^?4,  .00  00?!/ 

REAL* 9  NGI  5)/2,  ,  Ra  ,l«,  ,2  5,  fQ,/ 

RCAI  *8  TS{c)/,o?o?t,'>n44t.,  0007?,,on*062,,'^00m/ 

REAL*R  NSf  *)/?„  ,o,  ,1°,  ,!««,  ,4./ 

RFAL*R  IC(  ■>  )/,  oorqr  t  ,ooooot  ,  nnn^R/ 

REAL*R  NC(3)/2, ,8„ , 1», / 
RCAL*B  T A( ?)/,001 56, >OO0^RQ/ 
REALMS  NA(  ■>  )/?,,  «,  / 
NP0TL=6 
NPQTG=5 

npots=e 

NP0TC=3 

NPOTA=? 

SB  IS  THF  SUMMATION  OF  IONIZATION  POTENTIALS  F  0°  THE 
EXPERIMENTAL  ATRM, 

SR  =  0» 
XMASE  IS  THE  REST  MASS  OE  AN  ELECTRON  IN  mfv, 

XMASE=0  511006 
BETA  TS  THE  NORMAL  V/C  USED  IN  RELATIVISTTC  CALCULATIONS* 
GAMMA  IS  THE  NORMAL  TFRM  USED  I N'  RELATIVISTIC  CALCULATIONS 

BETA=DSORTUE!  '-'El  +  ^o^EP  XMASE)  /(EI  *EI+?,  *c  I*  XMASE 

1+XMASE  ^XMASE)  ) 
GAMMA=U/( l.-BETA**2)**(l*/2o) 

AT  THIS  °OINTt  thf  PRnr.RAM  MiJST  ehmphte  the  proper  CON- 
STANTS FOR  THF  FXPFPIMCNTAL  Z„   A  SFPIFS  ^e  LOGIC  STATF- 
MENTS  SELECT  THr  CORRECT  FORMULAE, 

13  IF(7«GT,13*l)GOT02c> 

COMPUTE  SB  FOR  cxpcR TMpNTAL  Z,   B  IS  THE  IONIZATION  POTEN- 
TIAL FROM  n^c  SHELL. 

DO  71  I=1,NPHTA 

B=IA(I)J'NA(nj'nLOG(">^cI/{IA(n-<'(l,-B,rTA*^?))) 
71  SB=SB+B 

ESTABLISH  VALUES     FOP     R?     AND    Al  ,        THESF    A p c    MATFPIAI      DE- 
PENDENT TONSTANTS    USrD    TO    COMPUTE     THE     AV^AQ^    FM«"pr,Y 
LHSC     BY  T0MI7ATinN    AMD    A1??     SUPPLTCD     BY    Tf|r     <;TFPNHEIMCp 
PAP£P,  THESE     ARE     B     AND     A,     R«=SPCCT  I V  El.  Y  ,     EOPM    HIC     PAPED, 

r?  =  ] ft,77 
Al =.07^0 

OFLT  IS  A  PART  OF  THE  EXPRESSION)  BY  cTpoNHcIMrp  cno  THE 


73 


c 
c 
r. 
c 

c 
c 
c 

c 

c 
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AVFPAGE    FNCRGY    !.nSS    DUC    TO    ionization,        NiJMRFRS    Hedptn 
PPPPFSFNT     OTHER     STFPMHPIMFO     CONSTANTS    AS     FOLLOWS: 

SMALL     A        „oqr>6 

SMALL    w       3*510 

X!  3.000 

-C  4*210 

0ELT=4-,  fi^  +  nLnr,i  oj  bcta*GAMma)-4,21  ♦-,  0906* 
1  (3„-nL0G10(RFTA"GAMMA)  )"*3 , ^1 
GO    TO    11! 
?Q     IF(7.GT0^q,! )GnTn^O 
on    72    I=1,Npotc 
B  =  IC  U  )*NC(n*DL0G(2*  *EI/{  ICC  I  )*( 1„ -BETA**? )  )  ) 

72  SR=SR+R 
B?=!  5,oq 
Al=.07^1 

DELT=4„606*OLOG10(  BETA*GAMMA  )-4-,74+a  11QQ* 
H3«-0LQG10<  PFTA*GAMMA)  )**3,^8 
GO    TO    ill 
50     IF(7eGT,50,1 )G0TO64 
00    7?    1=] ,NPOTS 
B=ISU  >*NS<  I  )*DLOG(  ?**-EI/<  ISin  +  lh-RFTA^?)  I  ) 

73  SB=SB+B 
B?  =  13»  93 
Al=,0647 

OFLT=4,606''OLOG10(  RET  A  *GAMMA  )  -6  -,  2R  +3  404*  (  3  ,  - 
10L0G10( BETA' GAMMA)  ) 

GO    TO    111 
64    IF(Z,GT,64^nGnTnfl? 
DO    74    I=l,NPOTG 
B=TG(I  )*NG< I >*DLOG(?,*Fl/(  IG(I)*(1 , -BFta**2)  )  ) 

74  SR=SP+R 
B2=13.3 
Al=,06?4 

OfLT=4-,  606*DLOG1  0(  R^TA*GAMM  a  )  -  6-,  8+<,418*<  3  .  - 
IOLOG10( BETAJGAMMA)  ) 
1+*?. 10 

GO    TO    111 
R2     TF(7.GT,R?„1  JGPTnqnrj 

DO    75    I="»  ,  NPOTL 

B=IL (I )*NL( T )*DLOG(?,*FI/(  IL ( I)  * (1  , -BETA**  2 ) )  ) 

75  SB=SB+B 
B  2=12.81 

Al=,06OR 

DFLT=49606*OLGG10( BETA^GAMMA  )-6  ,0  3  +  ,  0  65  2^- 
lf  4.-DL0G10C  RETA*GAMMA  J  )**3,41 
GO    TO    111 

AR     IS    TH^    QUANTITY    SMALI      A,     FROM    R|  UNCK    ANn    WcSTPHALt 
MULTIPLIED    BY    R,     the    TARGFT    THICKNFSS    IN    CENTIMETERS, 

111     AR=0, 154*7 *T/A/RFTA/RFTA 

COMPUTE    RS?,    TH^    SMALI      R    SQUARED    FROM    THE    BLUNFK    AND 
LFIS^GANG    PAPER, 

BS2=3.*SB/( 7*AR) 

VAP1,     VAR?,     ANO    VAR"*    A  R c    CONVENIENCE    STORAGE    LOCATIONS 

USFO    TO    STORE    COMPUTcO    PORTIONS    Hf"    S  T  rR  MHC  I  M  FR  '  S    4Vct?^c 
ENERGY    LOSS, 

VAR1=A1*T/BETA**2 

VAR7=R?  +  A4^  +  ?„*DL0G(RETAJ- GAMMA) 
VAR3=DL0G( CI )-RETA**? 

QAVF     IS     THF     FINAL    VALUE    FOR    THF     AVERAGE    FN^PGY    LOSS. 
Q^V^:  =  VAR1*(V^R'  +  VAR:,-DEL'^) 
B1=DSQRT(BS?J 

X     IS    THF    PORTION    OF     SOME    TOTAL     ENERGY    LOSS    Q    OIF     T^ 
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c 

c 
C 
c 
c 
c 
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c 
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c 

r 
c 
c 

r 

c 
c 
c 
c 
c 
c 

r 

f. 

c 
c 

r. 
c 
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PEAL  FUNCTION  WR  An  *<M  X  f  p  I ,  T  ,  7  t  A  ,  0  ) 

pljPPDSF:     THIS    FUNCTION     CALCULATES    THAT    PART     or    THF    cNFROY 
0  ISTRIPtJT  I  ON     nF     AN     F|  PfT^N     HUE     Tn    LOSS     IF     PW^r,Y    «Y 
RAOTATI0N    WHILF    PASSING    THROUGH    AN    ARSORRrp, 

IMPLICIT    RcAt  +  R(  A-H,n-7  ) 

U     IS    AN    TNTFPMEDIATF    STORAGE    LOCATION    USFO    IN    CnMPUTIMH 
ALPHA* 

U=183.0/Z**U. /390) 


COMdijtc    VALUE    or    AtRHA^R    fop    TAPGCT    THTCKNCSS,     ATOvjr 
NUMBFR     AMO    ATiVIC     WEIGHT    HP    ^XP^R  I MENT  Al      ars^RCr, 
THE    TARGFT    THICKNESS     FXPRESSFD    TN    CENT  I m^TFP s , 

ALPHR^^OO1  4*T*Z*7/A*{  4„0/3.>0*QL0G(U)+!  ,n/°>°) 

COMPUTE  VALUE  for  pf  a  NORMALIZING  CONSTANT, 

R  =  K  /DGAMMA(ALPHP+!»  ) 


IS 


OF 


COMPUTP  WRAO,  TH^  NIJMRFP  np  COUNTS  cXPPCTFn  AS  A  RFSUI 
PAOTATION  PNCPGY  LnSSESo   X  PcPRcSrNTS  THAT  pni<TT™N  0«" 
Spmt  TOTAL  ENERGY  LOSS  0  WHICH  IS  NOT  LHST  BY  IONIZATION 

WRAO=( R*ALPHPt (X**ALPHO )  / ( F I ** ALPHP )) /X 

RETURN 
ENO 


RFAL    FUNCTION    BFAM*8 ( FI ♦ HW, X ) 
IMPLICIT    PFAL*«( A-HfO-Z) 
Y=?, 

ALFA  =4,  *OLOGC  Y    )/(HWfcHW) 
BEA^=D'XP(-(  IHII*"I  ^ALFA) 
RETURN 

ENn 


RCAL    FUNCTION    FCT*  8  (  X»  7 1  A»  T,  E  \  ,  o  ) 

PURPOSE:  THTS  FUNCTION  mijl  t  i  pl  I  0  AT  T  VH  Y  Jhjms  TOGETHER  thc 
PRFnfCTinN  OF  c^cp^y  |O^S  DUF  TO  IONIZATION  AMO  TMP  LOSS 
DUE     TQ    RADIATION, 

I^PLirTT    RFAL*«( A-HtO-71 

FCT=WRAn(x,FI  ,T,  7,  A.Q     )  *  WI  ON  (  X  ,  F  T  ,  T  f  7  ,  A  ,  0     ) 
RFTU°N 

END 
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SlRRntJTTNP     GRAPH  (THX.THV,  FXX,CXY,T!  ,TM^,Ti,HWf  n  ,T, 
ICORRfCTSNORtQP) 

PUPPOSc:     THIS    SURPnUTINc    PL^TS    THF    RFSULTS    oc     TH^npy    CAL- 
CULATIONS    AND    FXP^PI MENTAL     RESULTS. 


IMPLICIT    rcai *Q( A-H,n-7) 
RFAl *P     TITL^Xt 1  )  /• mfv  «  / 

R^AL*3     THLrY(1)/»M1RMALI7Fn 
RC&I  ^q 


COUNTS 


i  i 


DIMENSION    T1(^),T?<4|,TM^),TM^I,THX(1^)  ,THY(30O) , 
1FXX(?00)  ,PXY{?oo|  ,  CTSUfM?^  )  »C  TSDNHOO ) 


TTTL^Ot  ?)  /»  Roy 


DFTERMINF  SCALING  FACTORS  AMD  CHANGE  DATA  Tn  DISPLACEMENT 

X=8, 

Y=5. 

SIXIN=(CI-QP)*?) 

ISIXIN  =  S!XIN  +  ,e; 

SIXIN?= I  SI XIN 

SIXIN=SIXIN?/2> 

IF(HW. EO.0JGOT090 

XNCR2>i 5*HW 

INCR2=XNCR?+*5 

IFCINC2»LT«1  )INCR?=1 

XINCR2=INCR2 

XINCR=XINCR7/2^ 

SP=SIXIN-6,*XINCR 

IF(SP«LT.O)GOT090 

G0T091 
oo  SP=n 

INCR  =  (EI-l<.5*T)/6.+o5 

XINCR=INCP 
91  N=0 

1=0 

20  1=1+1 

TF(EXX(  I  |,GT,1, )  GOT021 
GOTO?? 

21  N=N+1 
GOTO  ?n 

22  IF(N.F0,0)GOTO100 
DO  23  J  =  1,N 

EXX< J)=<EXX( J1-SP) /XINCR+CORR/XINCP 

IF(J.EQ.N)G0T073 

CONT INUE 


)  GOTO  31 


C 

c 


23 

100  1=0 
M  =  0 

30  1=1+1 
IFCTHX!  I)«GT.l 
GOT032 

31  M=M+1 
G0T030 

37  On  3^  J=lfM 

THX( J)=(THX(J)-SP)/XINCR 

IF(  j0ro»M)r.cTn^? 

IF(THY(  J+l  KGT  -,TUY(  J)  )  THYM  A  X  =  THY  (  J  +1  ) 

33  CONTINUE 

IF(N.EO.O) GHT010! 

SCALF=4„ 

IF(T.GT04„ R)  SCALF  =  3o 

on  40  j=i,n 

FXY<J1=EXY(J)/CTSN0R*SCALE 

^0  CONTINUF 

101  DO  50  J=^fM 
SCALE**.. 

IF(T^GT,4,  q)  SCA!.F  =  ^, 
THY( J)=THY(J)/THYMAX*SCALE 

^0  CONTINUE 

SORT  EXX  AND  ELIMINATE  VALUES  OFF  GRAPH 
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60 


62 

61 
10? 


7? 

71 


IF( 

K  =  N 

on 

I  PI 

on 

Ic( 

TCM 
TCM 
PXX 
CXY 
FXX 
EXY 
COM 
J  =  0 
KK  = 
K  =  N 

on 

TF< 
IF( 
J  =  J 

N  =  J 
EXX 
CXY 
1F( 
GOT 
N=J 
KK  = 
CON 
KK  = 
J=0 
K  =  M 

on 

IF( 
IF( 

J  =  J 
M  =  J 

THX 
THY 
IF( 
GOT 
M=J 
KK  = 
CON 


»N 

c,  PXX(J  )  ) 


J) 


NoEQ,^) GOTniO' 

-1 

60  1=1 ,k 

=  T  +  1 

6  0  J  = I P 1 

FXX( T  ),l 

P=FYX( T  ) 

=  PXY(  T  ) 

( I)=CXX( 

( I)=rXY( J) 

(  J)=Tf:MP 

i J)=T=M 

TINUE 


61  I=1,K 
KK,F0„1 IG0T061 
CXX(  I  )oLT„0)00T06i 
+  1 

( J)=^XX( I ) 

( J)=rXY( I ) 

CXX( J),GT„X)GOTnA? 

0  61 

1 
TINUE 

0 


71  I = 1 , K 

KK.FQ^l )0nTn71 
THX(  I  ),LT,0)r,nTn7i 
+  1 


GHTP  60 


( J)=THX( T I 
( J) -THY(T) 

THX(  I  ),,  GT, 

071 
-1 
1 
TINUC 


XIG0T07? 


DFTFRMINE  PRROR  BAR  MAGNITUDE 

TF(N,P0,0)GnTOin3 

OH  700  1=1 tN 

CTSON(  T)=FXY(I)-n<;0RT(rxY(l)A-?C^,)/?0o 

CTsUP(M  =  FXY(n  +  n<;QRT(FXYfn'-?Orv>)/?or 

7°0  CONTINUE 

INITIALIZE  "LOT  AND  WR  IT^  I  DENT  Tr  I  CAT  I  ">N 

103  CALL  PLOTS 

CALL  SYMBOL (0,  0,  ,?°,TTTLr^»r»l  6) 
CALL  PL0T(0,2.t-3) 

DRAW  OUTLINE  n^  GRAPH 


7  =  0 

CALL 

CALL 

CALL 

CALL 


PI  nT( 

PLOT( 
PLOT( 

PI  PT  ( 


»Y, 
»Y, 

,7, 
tZ, 


DRAW  ^UTLIMC  OF  TITLE  BOX 

Bl=„^ 
B?  =  ^  ,A 

84=4.6 
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70! 


5  =  R 
*  =  R 
ALL 
ALL 
ALL 
ALL 
ALL 
ALL 

5  =  R 
ALL 
^  =  R 

6  =  R 
XP  = 
ALL 

5  =  R 

6  =  P 
ALL 

5  =  R 

6  =  P 
ALL 
F(N) 
6=B 

5  =  R 

7  =  R 
ALL 
ALL 

6  =  B 

5  =  R 
A  =  B 
ALL 

6  =  R 
ALL 


!  +  >! 
3+,  I 

PI  0T{ 
PLHT( 
PI  ~>T( 
PI.  OT( 
PI  HT( 
SYMRH 

5+„  2^ 

Nll^RE 
54-1  „  32 
6  +  -,! 

NDMRP 

6-.1 

NUMRF 
5-2*76 
6+o  2^ 

SYMRO 
,EO*0) 
6+on7 
5+„  R2 
S+3  c  ^ 

PLOT 

PLHT 
6-. 07 
5-,  R? 
(S+o  ">4 

SYMRP 
6+,2^ 

SYMRO 


81 , B*,"M 
*?, RR, ">) 
R?,R4,  ?) 
B1,B4,?) 
R1 I,  R3,  2) 
L(B5,B6,,14,  Tlf0>^,3^) 

P(B«5,B6,„14,T,0,0,3) 


D (85,86, J07,FXP,0o0,-1 ) 


R(B5,R6,.,14,FT  ,0,0,2  ) 


L  (  R  5  ,  P  6  ,  ,  1 4 , 
GDT0701 


(85,86,3) 

(P7,RA,2) 


L(B5»B6,ol4, 
L  (B5,R6,,  1.4, 


PL^T    TIC    MARKS 


1  S7n 


l  5  BO 


lf>70 


1  690 


1K71 


!5Q1 


=  Y-,1 
ALL    PL 
ALL    PL 
=  G  +  ,1 
p    <G„L 
=  X-»1 

ALL    PL 
ALL    PL 
=  G-»l 
F(G.GT 

=  o2 

ALL  PL 
ALL  PL 
=  G+1. 
F  (  G  ,  L  T 
=Y-.2 
ALL  PL 
ALL  PL 
=  G-1„ 
F(G,GT 
=  •1 
=  Y-,1 
ALL  PL 
ALL  PI. 
=  H-»! 
F(H.GT 

=  ,1 

ALL  PL 
ALL  PL 
=  H+.l 
F(H„LT 
=  X-,2 
=  Y-1, 


0T( 


G  ,  Y  »  3  ) 
G,H, 2) 


F..X1G0  TO  1S70 


0T(G,H,3) 

0T(G,0o  t2) 

•0. )G0  TQ  1580 


0T(G,H,3) 
QT(G,0,2) 

,  X)  GO  T0  1670 


0T( 

0T( 


G,H,R) 
G,Y,?) 


Oo  )  GO  TO  16R0 


0T( 
0T( 


G,H, 3) 
0o,H,2) 


,0.  )  GO  TO  1 *71 


GT( 
1T( 


G ,  H  ♦  3  ) 

X,H, 3) 


T?  ,  0,0,3^) 


T3,0»0,32) 
T4, 0,0, 32) 


Y)  GO  Tn  15B1 
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167]  CALL  PlOT(G,Hf^l 
CALL  PLnT(X,H??) 

IF(HoGT.O)     OH     Tn    ]A7i 
G-m? 

H=l. 

!**!     CALL    PLnT(GfHt">) 
CALL    PLGT(0,Hf^>) 

IF(HaLT,Y)     GO    TT     l^qi 

c 

C    LAB^L    AXIS 

C 

CALL    SYMRPL(3,5,-^t,}'MTTTiexf%4) 
CALL    SYWRnL(-e?e;f1,'56>t>14,TTTLFY,o?,,?4» 

C 

C    NUMRFR    ENERGY    AXIS 

C 
1050    G=-e2P 

FLO=SP 
1^60    CALL    NMMRPPCG,  H,,l^fFLn,n,->) 

G=G+U 

IF(G0GT,X)     GH    Tn    1070 

FLO=FLO+XJNCP 

GO    TO    1^60 
C 
C     PLOT    FXPERIMFNTAL    OATA 

c 

1070   if(n,F0o0)  nnmioA 

CALL    LINEC EXX, Eyy, Nf?,-M 


C 

c 

C 


c 
c 
c 

c 
c 
c 


PLOT    ERROR     FLAG^ 

nn    101  o    1=1  f\i 
ES!=FXX( I >-*05 
c5?  =  cXx(  I  >  +  ,  r>  r 
CALL     PLOT(cSl,CTSnP(  I)  ,3  ) 
CALL    PI  nT(rS?,CTSUo(  I  )  t?) 
CALL    PLOT(FXX(n  ,CTSUP<  I  )  ,3  ) 
CALL    PLnT(EXX( I) ,CTSOM(  I  ),?) 
CALL    PL0T(FS1 ,CTSON( T) ,3) 
CALL    PLGT( ES?,CTSOM< I) ,3) 
1010    CTNTINUC 


PLOT    THEORY    CURVE 

104    CALL    LINE(THX, THY, M, 2,11 

FINALIZE  PLOT 

CALL    PLDT(0«,,8a,-3) 

CALL    PLOTF 

RFTIJRN 

END 
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